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^ : Abstract 

Resonance-mediated many body decays of heavy mesons are analyzed. We 
\ focus on some particular processes, in which the available phase space for 

ly-^ I the decay of the intermediate resonance is very narrow. It is shown that the 

'nI" ■ mass selection criteria used in several experimental studies of D and B meson 

decays could lead to a significant underestimation of branching ratios. 



^ Three or more body decays of heavy mesons are usually dominated by intermediate res- 

p • onances. The latter are short lived states that cannot be directly observed: only daughter 

p • particles produced through their decays reach the detectors The detected final state 
D ! results from the interference of all possible intermediate channels. Thus, one has to disen- 
^ ' tangle the quantum interference to understand the underlying physics. A powerful technique 

to split up the various resonant channels is the so-called Dalitz plot fit 
^ ■ In general, if a given resonance proceeds within a kinematic region where no other res- 

. onances occur, interference effects are assumed to be negligible. In this case, the branching 
ratio of the heavy meson decay to that particular resonance can be measured in a simple 
way: one just counts the number of events in which the invariant mass of the decay products 
lies within a small window around the resonance central mass. The size of this window is 
chosen according to the width of the resonance. This is a natural and thus widely used 
measurement technique 0. 

In this Letter we discuss the validity of this method. In fact, the resonance is a virtual 
state, and its squared four-momentum can reach in principle any value within the allowed 
phase space. This is certainly not a new statement. The interesting fact is that, as we will 
see, for some particular decays one could find a very large number of (naively, unexpected) 
events such that, even if they proceed through a particular resonance, the invariant mass of 
the detected particles is indeed very far from the resonance mass shell (with "far" , we mean in 
comparison with the resonance width). These events would be missed if the counting method 
simply considers a narrow window around the resonance central mass. As shown below, this 
is the case for some processes in which the resonance decay rate is kinematically suppressed. 
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We are aware that the physics involved in this discussion is well known. However, we 
believe that the magnitude of the effect has not been well appreciated so far, and the 
described counting method is not always safely applicable. Here we present simulations of 
actual decays and discuss the significance of the effect as well as the expected experimental 
difficulties. 

Let us first describe those aspects of many body decays relevant to our discussion. Con- 
sider a heavy meson P decaying into a final state given by three detected particles, di,d2 
and (is, and assume that the decay proceeds through intermediate resonances Ri,R2, etcQ. 
The total amplitude of the decay is the sum of the amplitudes of all partial channels, each 
one mediated by a resonance Rf 

Atot = ^^AR^ . (1) 

For simplicity, let us assume that the decay is dominated by a single resonance R, 
in such a way that A — Ar = A{P Rd^, R did2). It is natural to describe the 
process by considering three independent stages: resonance production P — >• Rd^, resonance 
propagation, and resonance decay R did2. The amplitude factorizes then as 

A{P Rd3 rfic/a^s) = 
A{P Rds) X BWr,i2 X A{R rfirfa) • (2) 

The amplitudes A{P — > Rd^) and A{R — > did2) have to take into account the con- 
servation of angular momentum in the decays, as well as the energy dependence, usually 
parameterized through form factors (we will come back to these two ingredients later). As 
usual, we describe the resonance propagation by means of a relativistic Breit-Wigner func- 
tion § 

BWR,i2iml,) = — 1 — - , (3) 

where mo is the resonance mass, F is the resonance width, and is the invariant mass of 
the outgoing particles di and 6/2, rn\2 = {pi + ^2)^- 

The function \BW{ml2)\ is peaked around the resonance mass, and decreases with a 
rate given by F. This behavior reflects the fact that i? is a virtual particle that — according 
to quantum mechanics — can have any invariant mass m^25 relative probability of each 
"mass" being weighted by the factor (777.^2) It is then natural to measure BR{P —>■ 
Rd^) by simply counting the number of detected particles di and d2 for which the value 



of \jm\2 liss within a region of the order of (mo — F,mo + F). This technique has in fact 
been used for many years 0. However, as we will show in the following, the usage of this 
approach is not always safe. 

The decay of a resonance is usually driven by the strong interaction. Thus, resonances 
have relatively large widths — some dozen or even some hundreds of MeVs [D]. However, 



The decay also proceeds through a non-resonant channel, which is not relevant to our discussion. 
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if a particular decay is somehow suppressed, the resonance may have a longer life, and its 
width can be as small as some MeVs or less. This happens in particular when the resonance 
central mass mo is very close to the threshold of its decay to did2- In this case, the phase 
space available for the decay turns out to be very narrow, and the (virtual) resonance could 
be allowed to decay through other channels which were in principle expected to be strongly 
suppressed in comparison with the "natural" strong channel R — > did2- 

This is the case, for instance, of the 0(1020) vector meson. For this resonance, the natural 
decay channel is KK, in the same way as the natural decay for p is two pions. Nevertheless, 
due to the small phase space available — 32 MeV and 24 MeV for the charged and neutral 
kaons, respectively — the corresponding (p branching ratios are "only" 49% and 34% for 
K^K~^ and KK, respectively. Electromagnetic decays, that have branching fractions as 
small as 10~^ in the p decay pattern, are of the order of 1% in the case. Accordingly, the 
width is about 35 times smaller than the p one. 

Other examples are low mass D* vector mesons. Their natural decay channel is Dvr, in the 
same way as the natural decay for K* is i^vr. But here, the phase space is as small as 7 MeV 
for the L)*°(2007) D^n^, 6 MeV for both D*+(2010) D°n+ and L)*+(2010) D+tx^, 
and 8 MeV for D*^ — > D^tt^. As a consequence, measured branching ratios of these decays 
— which otherwise should reach almost 100% — are as "small" as 62%, 68%, 31% and 6%, 
respectively^]. In the D*°(2007) decay pattern, the electromagnetic decay is as large as 
38%, i.e., of the same order of magnitude of the strong one (in contrast, in the K* case, 
electromagnetic decays are of the order of 10~^). Accordingly, resonance widths are quite 
small: the width of D*~^{2010) has recently been reported to be as small as 0.1 MeV 0, 
whereas for D*°(2007) and D*^ only upper limits are known, presently of the order of 2 
MeV. 

Let us face the study of heavy meson decays mediated by these particular spin one 
resonances, focusing on cases in which the detected final state includes their natural, and 
yet highly suppressed, strong channels. We describe here a usual situation [|I|, where both 
the initial heavy meson and the particles in the final state are scalars, P = D,B, and 
di = n, K, D. In this case, Eq. (H) can be conveniently written as 

AiP -y Rd3 rfirfarfs) = Fp^Rd^ Fn^a^d^ (-2pi ■ Ps) BWr^i2 , (4) 

where Fp^Rd^ and FFt,did2 are form factors, and the three-momenta pi and ps are evaluated 
in the resonance rest frame. The explicit momentum dependence in (^ follows from Eq. 
(0), just assuming Lorentz invariance and summing over all possible polarizations of the 
intermediate vector meson resonance. The differential decay width of this reaction can be 
written as 

^ (2^32^''^'''^'"'2 "^""'^^ ' 
where m|- = {pi + Pj)^ and M is the mass of the decaying meson P. 



"^Df* Dfi:^ is an isospin violating decay (see Ref. j^). It thus has another, phase space 
independent suppression. 
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For the decays considered here, there is a strong suppression at the resonant peak, i.e. 
when — ""^o- This suppression is due to purely kinematic effects. To see this, let us 
take A ~ constant for a given value of — which means to neglect the dynamics of the 
decay — and integrate over the variable mf^ within the kinematic limits of the three body 
phase space. It is easy to see that 

^oci^nfti, (6) 

where pi is the three-momentum of di in the resonance rest frame. Since we are assuming 
that the mass of the resonance is just above the threshold, mo — mi + m2, at the resonance 
peak both particles di and d2 will be produced almost at rest. Thus Eq. (|^) implies a 
suppression in the partial width. 

The effect is even stronger in the particular case of a vector resonance. In that case, as 
stated in Eq. (|^), the amplitude is proportional to pi-p^. Assuming that the form factors are 
slow-varying functions of the phase space variables, the partial width is finally expected to 
be suppressed by a factor (|pi|/A)^, where A is some natural scale of the process, typically 
of order M. 

Eqs. (^) and (|^), which are certainly very well known, show up the main point of this 
Letter. The decay width T{P — >• ^1^2(^3) is driven by two competitive effects. On the one 
hand, the BW propagator strongly enhances the decay amplitude in the vicinity of the 
resonance mass. On the other hand, kinematic effects suppress the differential width at the 
resonance peak, hence the decay rate for other kinematic regions is comparatively enhanced. 
The usual suppression of the differential decay width for P — > Rd^ did2d^ outside the 
window allowed by the BW function is not obvious in this case. 

To clarify our point, let us emphasize that we are not claiming that there is an en- 
hancement of the resonant production probability P —>■ Rds outside the BW peak. On the 
contrary, the point is that due to the suppression of the resonance decay rate R —>■ did2 at 
the resonance central mass, the combined production + decay probabihty within and outside 
the peak could be of comparable orders. In other words, the total width (i.e., integrated 
over the whole phase space) is indeed small; the important point is that, even if the decay 
proceeds through a BW-described resonance, the relative weight for the decay rate near or 
far from the resonance mass shell is not just driven by the BW function. 

In order to show that this effect is not a simple academic thought, we present in the 
following a simulation of an actual process. We first present the pure theoretical estimate; 
afterwards, we will consider the experimental difficulties. Let us consider the decay 
£)*od+- £)*o _^ D^Tc^. Using Eq. (Q), it is possible to get an estimate for the differential 
decay rate dT as a function of m^g- Since the form factors are usually smooth functions, we 
will assume as a first guess that they are constant (experimental analyses show that form 
factor shapes have no significant effect on the total systematic error of a given Dalitz plot 
fit 0). In this way we can calculate the ratio 

\A\^d<^ 
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where (i$ is an element of the three body phase space, mo is the D*^ resonance mass 
(mo = 2007 MeV), F is the D*^ width, and n is a real number. F is presently unknown, its 
upper limit being 2.1 MeV with a 90% confidence level |^. The integral in the denominator 
is calculated over the whole phase space, while that in the numerator is limited to a window 
in m^2- Thus, r is a measure of the relative number of events that are expected to fall within 
the resonance peak. 

We quote in Table I the values of r for some input values of n and F. Our results show 
that the effect we are describing can be very strong if the resonance width is as large as 2 
MeV, and it remains quite significant even if F is of the order of 0.1 MeV. For comparison, 
we include in the last column the values of r corresponding to a fictitious resonance having 
the same quantum numbers as -D*° but a higher mass, mo = 2.6 GeV. This particle would 
not suffer the kinematic suppression in its decay to D^n^, consequently a small value of n is 
enough to get r above 90%. We have found that in this case the results for r are independent 
of the specific value of F. 

Figures la and lb show the kinematic distribution of the events for the three body 
decay —>■ D^ir^Df, assuming that the process is dominated by the Z)*°(2007) resonance 
channel^ and considering a D*^ width F = 1 MeV. The plots correspond to a Monte Carlo 
simulation of 10000 events, performed using Eqs. (^) and (^). Figure la is the Dalitz plot 
of the decay as a function of the invariant masses m^o^o and m^+^g, while in Figure lb we 
represent an histogram of the number of events as a function of m|,o^o- 

It is seen that the Dalitz plot shape in Fig. la is quite different from that naively expected 
for a decay mediated by a F = 1 MeV vector resonance. This becomes evident by looking at 
Figure Ic, where we show a Monte Carlo simulation of the same process, now shifting the 
mass to the fictitious value of 2.6 GeV considered in Table I. The striking difference between 
the event distribution in both plots arises from the kinematic suppression discussed above. 
However, the situation displayed in Fig. la can be misleading if one just looks at the events 
for which m^o^o is near the D*^ mass. In fact, despite the spreading of events along the 
whole phase space, the event density is still much larger in the peak region than anywhere 
else. Therefore, Figure lb could be mistaken for a F = 1 MeV Breit-Wigner function, with 
some background in the right sideband p. It would be then natural to apply the usual 
method, that means to consider that the events within a window of a few F around the peak 
amount almost the total number of decay events. But, according to our simulation, this 
would be wrong by a factor as large as 3 to 4 (see Table I). Indeed, even if the amplitude is 
peaked around the resonance mass, the phase space area outside the peak is comparatively 
so large that the number of events falling in this region becomes very important. 

Let us now discuss the experimental difficulties that could appear when data from real 
events are analyzed. First of all, resolution in these experiments is usually larger than the 
width of the intermediate resonances involved . As a consequence, it would be impossible 
to access to a direct measurement of the ratio r, at least for small values of n. Accordingly, 



^In fact, the decay can also proceed through an intermediate resonance D*~^ , as well as other higher 
resonant states. The inclusion of these contributions in our analysis does not affect significantly 
the results. 
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the experimental histogram of Fig. lb would be wider, arising from the convolution of 
a narrow Breit-Wigner function (physical) and a Gaussian function with a larger width 
(resolution). Nevertheless, the bulk of our reasoning remains valid: a large amount of 
physical events may fall outside the peak. A second comment refers to the difficulty of 
developing a more careful analysis of data, in order to include the contribution of these 
events. Indeed, according to Table I, the amount of events could be quite important, but 
at the same time they appear to be spread in a very large region of the phase space. Then, 
only a relatively low number of events per bin would be found outside the peak, and they 
could be hardly disentangled from the background, no matter which is the function used to 
perform the corresponding fit. 

It is important to stress that our theoretical estimates rely on two basic assumptions. 
First, we have taken the form factors to be approximately constant along the phase space. 
Second, we have assumed that the Breit-Wigner shape remains valid far beyond a small 
region around the peak. Whereas the first hypothesis is quite natural and does not have a 
significant effect on our results, the second assumption can be hardly supported from the 
theoretical point of view. Our simulations are in this sense strongly model dependent, and 
this has to be kept in mind when looking at the numerical values presented in Table I. 

In any case, our simulations can be taken as a severe warning, indicating that many re- 
ported results can be spoiled by the effect described in this Letter. For some particular pro- 
cesses, the magnitude of the corrections could be quite significant, and this should be taken 
into account — at least — in the corresponding systematic errors. This includes the analysis 
of the decays B D*^d^, B D*^d^ and B D*+d3, with 4 = 7r°, r/, K, etc. 
These channels have been measured using D^ir^d^ — respectively D^ir^d^ and D^ir^d^ — 
as final states, and imposing a cut in the invariant mass ^|)0^o — respectively ^|)+^o and 
m^+^Q. According to the simple analysis presented here, in all these cases the effect would 

be of the order predicted in Table I. In particular, in the case of the D*^ resonance, whose 
width has recently been reported to be 0.1 MeV, the effect is expected to be of the order of 
30%. 

Finally, let us mention that meson production and decay measurements could also be 
affected by the effect described above. Since the resonance can be produced in D meson 
decays, a significant amount of data is presently available, and many Dalitz plot analyses 
have already been done. However, notice that the threshold of the decay — * KK is about 
25 MeV away from the mass, so that the available phase space is not as narrow as in 
the D* case. Performing a Monte Carlo simulation for the decay D ^ (pir similar to that 
described above for the process B D*Ds, we have found that the expected effect could 
be in this case as large as 30%. 
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TABLES 



n 


0.1 MeV 


0.5 MeV 


1 MeV 


2 MeV 


fictitious 


1 


61% 


30% 


18% 


11% 


70% 


3 


69% 


36% 


24% 


17% 


90% 


10 


73% 


41% 


30% 


25% 


97% 


30 


75% 


46% 


37% 


36% 


99.5% 



TABLE I. The ratio r as explained in the text, for different values of n (1 to 30) and F (0.1 to 
2 MeV). The last column shows the values of r for a fictitious resonance with mass 2.6 GeV — see 
text. 
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FIG. 1. Simulation of the decay B+ D*°Df; D*^ D^tt^, with T = 1 MeV and 10000 
generated events, (a) is the Dalitz plot in the plane [(p^o +P7ro)^, iPj:,+ +P-n-o)'^] ) whereas (b) is a 
projection on the {pQo +Pt:o)^ axis, (c) is same as (a), with the D*'^ mass shifted to 2.6 GeV 
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